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Abstract 

Charged pions produced in very high energy hadronic interactions might be the dominant source of 
cosmic neutrinos in the GeV-TeV range. Spectral energy power of n + radiation by high energy protons 
moving in strong magnetic fields typical for magnetars is determined with a semiclassical treatment of the 
£ — . effective pion-nucleon model. The main characteristics emerging from a saddle point approximation to the 

summation over the allowed range of Landau levels is a sharp lower cut: > 0.25E P . The magnitude of 
the spectral power agrees in this region with the synchrotron radiation spectra of neutral pions. 

(N . 

1 Charged pions as cosmic sources of high energy neutrinos 

Identification of the cosmic sources of (ultra)high energy neutrinos will make a big leap forward with the 
completion of the construction of ICECUBE [Tj, KM3NeT[2| and other observatories specifically designed for 
i— 1 1 neutrino observations, and also with the cosmic ray observatory P. Auger and its northern counterpart. At 
i-C ■ present theoretical arguments are put forward for neutrino production in gamma ray bursts 0, magnetars jH[S], 
blazars[B], pulsars[7] and compact binary objects[5]. In general objects where protons might be accelerated to 
Q ■ high enough energy to reach the threshold of pion production in various elementary reactions are worth to be 



o 
o 



X 



investigated. High energy neutrinos are expected to arise from pion and subsequent muon decays. 



The typical reaction chain in presence of hard enough 7 quanta is[9j 

p + 7 — * A + — > n + 7r + — * n + fi + + — ► ri + + v e + e + + D^. (1) 

Energetic 7 quanta might arise from synchrotron radiation of fast protons moving on circular orbit in strong 
magnetic field of magnetars, for instance. In the estimation of the neutrino luminosity possible effects on proton 
acceleration arising from the presence of strong magnetic field were taken into account. Another mechanism pro- 
posed for 7r + production involves high energy proton-neutron collisions, assuming that also the latter component 
is present in the environment of a gamma ray burst [3] 



p + n^n + n + n + — > 2n + fi + + — ► 2n + e + + v e + + . (2) 

The problem of neutral pion production as a result of synchrotron radiation of pion field strongly coupled to 
the protons in presence of strong magnetic field B: 



C|' p + B^p + TT° + B (3) 

a ; 

was investigated first forty years ago by Ginzburg and Zharkov [10] and reviewed by Ginzburg and Sirovatski[llj. 
The treatment followed closely the classical theory of synchrotron radiation. It was shown to reproduce the 
results of first order quantum perturbation theory applied to an effective nucleon-pion interaction Hamiltonian. 
A perturbative estimate for the p — > n+ir + transition probability under the influence of a strong electromagnetic 
wave (i.e. charged particles considered as Volkov states) was given at the same time by Zharkov [12]. 

At present much higher cosmic magnetic fields are known, in particular in mag netars, where B = 10 14 " 16 G 
was observed|13]. Substantial progress has been achieved also in the understanding of acceleration mechanisms. 
Therefore it is timely to evaluate the old formulae with the actual parameters. This task was partially done by 
Tokuhisa and KajinofTl]. but no reassessment was attempted to date for the case of charged pions. 
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The subject of this note is to present a semiclassical theory for charged pion radiation in strong magnetic 
field: 

p + B -> n + 7T+ + B (4) 

and estimate the luminosity of this radiation for the highest possible energy range. Since any individual proton 
contributes only a single positively charged pion into the radiation, it is hard to associate a classical charged pion 
source with a single proton. If one assumes a very intense proton beam, however, then the loss of a restricted 
number of protons should not mean a strong back reaction, e.g. the source remains approximately intact. This 
kind of approximation represents the main hypothesis of this work, which assumes unchanged proton source 
during the building up of a classical ir + radiation field around the source. 



2 The model 

The Lagrangian of charged pions coupled to a nuclear source density is the following: 

C = (dfj, - ieA ll )n + (d ll + ieA^)^ - m 2 n + i:~ - g^/2(p~ N ir + + p%n~) 



(5) 



where the nuclear source terms are p^ — ipj^n, pjj — in^p, and g 2 ~ 14.6. (Throughout this paper the unit 
system h = c = 1 is used, SI values for some important dimensional combinations will be explicitly given below). 
The Euler-Lagrange equations of the pion fields, in Coulomb gauge (Aq = 0, diAi = 0) read as 



[d 2 - d l d l ± 2 i eAidi + e 2 + m 2 ] tt* 



lV2p% 



(6) 



In presence of a static magnetic field along the V direction (B = (0, 0, B)), it is useful to introduce cylindrical 
coordinates. Then ^ takes the following form: 



1 d 



d \ Id 2 



1 



* pdp\dpj p 2 dej) 2 dz 2 d<j) 4 



-9 V2p± . 



(7) 



The particular retarded solution describing the radiation of charged pions for the above linear partial differential 
equations is constructed as 



^(p, 0, z, t) = -gV2 J dp'dcp l dz l dt l G ± {p, 0, z, t; p', 0', z', t')p%(p'<p', z', t'), 



(8) 



where the Green's functions G ± are determined by the eigenvalues and eigenfunctions of the differential operator 
of the left hand side of ([7|), 



C r[p, (j>,z,t;p ,<p ,z ,t) = 2^ 1 = {G ) . 



Writing the eigenfunctions in the form 



eB 



$ = exp [i (fup + kz — vt)] R(x) , x = ~~^P 2 j 



(9) 



(10) 



the corresponding eigenvalue problem is reduced to the eigenvalue problem of an ordinary differential operator, 



xR" + R' 



x p, p, 1 
4 + Ax ~ 2 ~ 2^B 



lo 2 - k 2 - m 2 ) 



R( x ) = -\R(x). 



(11) 



This is the well-known problem of the relativistic Landau-levels, with the normalized solution (see for example 
in POd] 22.6.16 and 22.2.12) 



^ \ / r/^iVL e- x/2 x W2 L l £ l («), A = eB(2n + \p\ -p + l)-oj 2 + k 2 + 

V r ( ?1 + ImI + 1) 



(12) 



where L„ (x) denotes the generalized Laguerre polynomial, and R£(x) form a complete orthogonal set for 
ormalized eigenfu; 



n = 0, 1, 2, 3, . . . . The normalized eigenfunctions (llOj) are written in the following explicit form 

'eB 



' eB 



exp [i (fi<t> + kz- cot)} R£ ( '—p 2 



(13) 
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and the corresponding eigenvalues are listed in (|12p. From (0) the Green's function is found after performing 
the oj integration: 



G+{ P Ar^t- P 'A',z',t') 



eB 

W) 



n— fi— — oC' 

l 



.4 



sin L4 n , M (t - *0] - *0> 



(14) 



where ^n, M ,fe = fc 



eB(2n+ -/x + 1) > 0. 



The source density is proportional to the transition matrix element from a Landau state of a proton to a 
plane wave neutron. We assume that a macroscopic circular proton current is present and neglect the decrease 
of its intensity due to the n + radiation. The corresponding classical proton density is given by 



p+(p, 0, z, t) = yi - ^ i s(j> _ R)5{z) £ i-e 1 " 



(15) 



where v p is the velocity of the proton (1 — v% = M 2 /E 2 = I/7 2 ) on the z = plane (we do not consider the 
motion of the proton along the z axis) and R, w a are the relativistic cyclotron radius and frequency of the 
proton. Substituting this source into ([8]) the charged meson radiation field is determined after performing the 
^integration: 



n + (p,(t>,z,t) 



eB ^— \ — isign (u) .. , . . . . 

■ 2^ —^f> ex P ~ + si s n {^z)yi n ^z)\ 



7 2ir 
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eBp 2 



(16) 



Here f2 n . M = y/ p 2 w 2 — m 2 — eB(2n + \p\ — p + 1). The allowed range of the summation is restricted by the 
positivity requirement of the argument of the root. The case of positive and negative values of p represent two 
distinct series in the sum = ^ + + 1 where the respective range of summation is given as follows: 



E 



V 'eB + m 2 < pu> s 
< n < 



2eB 



E 



yJ{eB/uj s ) 2 + eB + m 2 + eB/uj s < \p\uj s 

ft <: r, <: (M" 3 ) 2 -2e-B|M|-™ 2 -e-B _ 

U — " — 2eB ~ n - 



(17) 



The pion field is localized on Landau levels radially, therefore the energy flux is nonzero only along the z axis, 
e.g. one has to consider only the 'Oz' component of the energy momentum tensor, 



Tq z 



dir + ( dir 



Of 



\W) ~~dz~ 



(18) 



and the power of the radiation is given by computing the surface integrals of this flux on some distant ±z$ 
planes: 



= / / d(f>(T z(z = +z Q ) - T Qz (z = -zq)) = E 



eBRf 



2\ 1 2 



(19) 



3 The spectral power of the pion radiation 

The classical-quantum correspondence tells that the energy of the emitted pions is given by |/i|w a . Therefore 
the /i-sum should be converted into iS^dntegral if we wish to derive a formula for the spectral power. For this 
conversion one has to set a stage for the actual order of magnitude of the characteristic quantities. 
The frequency and the radius of the cyclotron trajectory of protons with energy E p are 



1 , 

hp "f z Bq oj s eB V e l B A 

where 7 = E p /M p and B Q = M 2 /e (in SI units M 2 c 3 /(eh) w 1.5-10 20 G). In $B5§ R% depends on the combination 

_ eBR 2 _ 1 2 -1 B 
~ 2 ~~ 2 B~ 
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The actual value of £ for a typical cosmic setting can be estimated by taking the values (in SI units) B = 10 15 G 
and 7 = 10 5 , which yield £ = 7.4 • 10 14 . 

Going from the /i-sum to the ^-integration the pion radiation power can be written in the following form 
(again in h = c = 1 units): 



9"eB 



/ d ^E' 



;(i^ /w -(0)' 



V 



■s/El - eB(2n + 1) - ?/i 2 

/ d£,rV-= V 7 (22) 

7e_ V _ eB (2n + 2E„/w s + 1) - m 2 



where from l[T7|l and |2U| 



e£\ 2 „ „ eB 



= VeB + rnP, £L = \ ( — +eB + m 2 + — > 2£ p . (23) 



The second integral is irrelevant due to £L > 2_E p (in the classical approximation without back reaction the 
limitation on the pion energy E„ < E p is imposed by hand!). The differential power of the radiation per unit 
pion energy is given by 

dE n ir^ojs ^ \JE1 -eB{2n+l) - m 2 
Since £ is a dimensionless (large) parameter, it is useful to express all quantities in proportion to £, 

~ - = "{I" ? i - 2/3^-£ ^ 2 P at (26) 

With these notations, from (fT2")) and (TT6")) . 

( fl »-""<«) 2 ° f(^Ti) ^-« W (i|{° { (0) 2 . (27) 

Since the change in A for one step in n is very small, can be written as £ J dA, and the following final 
expression is obtained: 

^r = 9 -^ / A+ dA - t F^,A,Q, (28) 

where 

F*(AA,0 . (<-'».(0) 2 . A + -2^i "' 2/M ;: B/J; ° ) ■ (29) 

The final task is to evaluate this integral. The difficulty is that the indices of the generalized Laguerre poly- 
nomials both are proportional to its argument £, which takes asymptotically large value. It appears somewhat 
discouraging that the integrand contains the very small factor e~^. Therefore the integration should be per- 
formed with great precaution. A saddle point representation will be described in the next subsection for the 
asymptotic regime of the Laguerre polynomials and shown to cancel exactly the suspicious small factor. 

4 Saddle point evaluation of the spectral tt + radiation power 

The key step in the evaluation of the spectral power is to find the appropriate asymptotic representation of 
L^i (£) and with this for the factor F({3, A,£) in the integrand of 
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Figure 1: On the left: For values of b — 2{3a and S — 2/?aA in the shaded region (I), the saddle points have 
non-zero real parts. Otherwise (in regions II and III) the saddle points are on the imaginary axis. The boundary 
of the shaded area on the right hand side (b = 25) corresponds to (|39p. On the right: The complex saddle points 
(0, z a , Zb) of q a and qb when (6, <5)G I, the original contour along the positive real axis (C) and the deformed 
integration contours (C a , Cb). 



The generalized Laguerre polynomials can be expressed with help of an integral of the Bessel functions (see 
22.10.14 in [T5]), 

L n(0 = ¥r —^e /2 e^ dfce- fc fc"+^ 2 J p (2 v ^|), Z = f- R - (30) 
T(n + 1) J 2a[3 

Using a convenient asymptotic form of the Bessel functions, F(P, A, £) can be written in the following 

form 



'Jo dz ^( ( 



,M<3a(z) I e (i?l(2) 



(31) 



where the expressions of q a ,b(z) and N' can be found in the Appendix. One finds that both terms of the 
integrand can be analytically continued into the complex z-plane. One can find its complex saddle points and 
can deform the integration path from the real axis into a curve of steepest descent in the complex z-plane which 
connects the real lower and the upper limiting points with the saddles (see the right hand side of figure [1]). The 
integral for large fi(= 2/?a£) is dominated by the result of the Gaussian integration around the tallest saddle. 
The location of the saddle points can be given in different ways for different regions of the (/?, A) parameter 
space (see the left hand side of figure [T]) . From the detailed saddle point analysis it can be concluded that the 
exponentially small N' factor in (f5Tj) is cancelled by the contribution of the saddle point only when (/?, A) are 
in the region I. The interested reader can reconstruct the steps of this argument by the details given in the 
Appendix. The result of the saddle point evaluation is the following: 



n 1/4 



4A - (2a/? - l) 2 



sin(^ 1 (/3,A) + $ 2 (/3,A) + 7r/4), 



(32) 



where 



2Aarctan 



^4A - (2(3a - 1) 
(2a/?- 1) 



^4A - (2a/? - 1) 



(4a/? + 2A)arctan 



^4A - (2a/? - 1) 
(2a/? + 1) 



J4A - (2a/? - l) 2 n JAA - (2a/? - l) 2 

$ 2 = arctan-^ \ _ 2arctan- v v 



(2a/?- 1) 



(2a/? + 1) 



(33) 
(34) 



when (/?, A)e I, and F(/3, A, £) rj when (/?, A)e II or III for large ^. 

Returning to the expression of the differential power of the radiation per unit pion energy, the factor 
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Figure 2: The differential radiation power per unit 7r + energy for different proton energies in comparison with 
the 7Tq radiation given by [14] for strong (B = 10 16 G) magnetic fields. 



F 2 (0, A, £) appearing in (|28p is just the square of the expression which appears in (f3"2")l . therefore 
dP g 2 M p £ f dA 



dE v IT 7 J ~ IJ 2 aA m*/M$+B/B 

X 4A - (2^ - 17 [1 + Sln (2e<&1 (/?Q ' A) + 2$2 (/3a ' A))] (35) 

The strongly oscillating term of the integrand would give 0(l/£) contribution therefore one can neglect it and 
one finds a simple result which is different from zero for a restricted range of E w /E p (0.25 < E^/Ep < 7 2 , for 
7 2 » 1) : 

dP _ g 2 M p f A+ dA 1 



dE. 



7 7(2/3a-l)V4 Jp2 aA ml/ml+B/B ,/ 4 A - (20a - 1) 

g 2 E n l b, 

rr, II — < — 

2tt 7 2 ' 4 ~ J5 



, ^ 7<^<7 2 - (36) 



This simple result can be used actually only for E v < E p since the energy of the radiated pion cannot exceed 
the energy of its parent proton. The result is compared next with the results of [2] and [TU] . 



5 Numerical results 

In figure [2] we plotted the above expression and the corresponding irg spectra for different proton energies at 
B = 10 16 G. In this plot we restrict the emitted pion energy to the physically interesting region (E„ < E p ) as in 
[T4] . It is interesting that the width of spectrum in our case is independent of the proton energy (in logarithmic 
scale) in contrast to the case of tto emission. This behaviour is obvious from the very simple expression (|36|) . 
and this is the consequence of the large £ asymptotic form of the Laguerre polynomial (|32|) . which becomes 
exponentially small when = E n /E p goes over from region I to region II or III (see figure [T]). The location of 



G 



the intersection of the boundary of region I and of the integration range of (|36|) is independent of £ (and B) and 
lies at E^/Ep = 1/4. In our approximation the slope of the spectrum variation is infinite at E^/Ep = 1/4 due 
to the very different behaviour of the Laguerre polynomials in regions I and III. As a matter of fact, our saddle 
point analysis is not valid in a very small vicinity of the boundary of regions I and II. due to the fact that the 
second derivative of q a ^ in ([31]) tends to zero when /3, A reach the boundary. The detailed investigation of the 
saddle point structure of Laguerre polynomials shows that the width of this "transition" region is ~ l/£, and 
the slope of the spectrum at E^/Ep = 1/4 in figured] is proportional to £ (logarithmic scale!). 

Another interesting property of the n + radiation is that the total emitted pion energy (integral of the 
expression (|36|) over E n ) is independent of the magnetic field in our approximation. This B independence is 
also shown in [T3] for the ir° radiation at large magnetic fields and/or proton energy, in contrast with [ID] , 
where the total emitted 7r° energy is proportional to B 2 . The origin of the difference between these two results 
(although both papers are based on the same semiclassical treatment of the ttq radiation, calculated in [ID]) is 
that in [2] the differential power of radiation is integrated over E v from to E p , while in [TDJ the integration 
range is infinite, and the dominant contribution of emitted energy comes from the unphysical (E„ > E p ) region 
of the spectrum in case of large B and/or E p . 



6 Discussion 

The prospective astrophysical interest of the very high energy neutrino source discussed in this paper depends on 
the simultanous occurrence of very strong magnetic fields and very energetic (e.g. 100 TeV) protons. Different 
propositions put proton acceleration at different distances from the surface of the neutron core, therefore the 
actual magnitude of the magnetic field strength might vary strongly. Most optimistic is the so-called "light- 
cylinder" accelerator mechanism, which in fast rotating objects can be located close to the surface. In this case 
the magnetic field strength is not much weaker than that directly measurable on the surface. Since surface 
keV-photons are necessarily present a concurrence between resonant photoproduction of charged pions (from 
the decay of A + ) and non-resonant synchrotron radiation- like-production Q should be envisaged. There could 
be 10-20 magnetars in our galaxy, with the great observational advantage that their activity is not restricted to 
only a few minutes, like in case of GRB's. Assuming some well-defined galactic distribution for these objects, one 
combines it with the neutrino luminosity of single objects to estimate the muon flux generated by ^-detection 
in ice. 
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Appendix 

We present below a detailed derivation of the asymptotic behaviour of ([3l]) when £ >> 1. The notation looks 
somewhat simpler if one goes back to the original Landau level wave functions 

F(P,A,0 = K(^), b^af3. (37) 

Relations between the variables A = n/^,fi — 2bt; were used in the main text. The integral representation of 
this function has now the following form: 



n+n/2+1 



K ( 2b) ^/T(n+l)T(n + f ,+ l) ^ ( 46 ) ( \ , 

x / dxexp(fj,[-xb/2+ (5+ 1/2) In x]) J^fiy/x) , (38) 
Jo 

where a new integration variable is introduced instead of the variable appearing in (|30|) and also a new parameter 
is defined: 

k = x&, Q< 8^< h -. (39) 
2 (i 2 y ' 
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In our application /i is large, one can use the asymptotic form of the Bessel functions given by 9.3.2 and 9.3.3 

of ma, 



1 



'2nfi 



7! J= cos(At(tan^i - K t ) - tt/4), if x > 1 
^tanh K2 exp(/i(tanhi^ 2 - K 2 )), if x < 1 



where 



cosisfi = y/l/x, 



cosh if 2 



175, 



tan if i = \/x — 1, tanhA^ = \ x , 

Therefore (|38p is broken up into the sum of two integrals: 



where 



h 



oo 2 

do; 



, Jx^vp 1 exp v 



^^J=iV(/ 1+ / 2 ), 



-— + (5+-) In x 



(40) 
(41) 

(42) 



X COS [ fl 

»1 



da; 



vi- T — acos 
1 



1 



AT 



1 



7T 

47' 

1 , 



~ + (S + ^) lnx + Vl — a; — acosh ( — ^= 



y/2-KiiT{n+ l)T{n + fi + l) 
Moreover, N can be replaced in the asymptotic regime using Stirling's formula for the T-functions, 



(43) 
(44) 
(45) 



N 



1 



n In 



exp 



1 1 

46 + 2 

fib 



2y/n + ly/n + /i + 1 



n + 1 + | In 



In 



2(ra + ^+ 1) 
fib 



2(n + fi + l) 1,A (n + l) lf \ 



(46) 



Now we can proceed with the evaluation of I\. The method will be similar also for I2, and will be just stated 
after the discussion of the case of I\ . One changes once more the integration variable in 



/, = 2. xi. ( -y ) / <1:\ : 



where one introduces 
and 



exp (^fiq a (z) - i~) + exp (jxqb(z) + i~ 
1 (a;- I) 1 / 4 = 



3a (*0 
9i>(z) 



z - + [8 + - ln(z 2 + 1) + Wz 2 - i acos 



z 2 + 1' 



+ ^ + 2") ln( ~ 2:2 + ^ ~ i ^ + 1 aC ° S \/ Z^+T' 



(47) 
(48) 

(49) 
(50) 



In (|47p the integrand can be continued analytically into the complex z plane, which allows the use of the method 
of steepest descent and stationary phase for the evaluation of the integral (saddle point approximation in the 
complex plane). For Rez > 0, 



q a ,b{z) = -z A - + [5 + - ln(z" + 1) ± iz T - In 



1 , 1 + iz 



2 1 - iz 



The saddle points are determined by the zeros of the derivative q' a b (z) 



Qa.b( z ) = ~ bz + 



(2(5 + 1)0 
z 2 + l 



±i T i 



z 2 + l 



(51) 



(52) 
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(6, 5)ei (M)en (6, S) e ill 

z a : {0, Zl } {0,z 2 }, (\z 2 \<l) {0,z 2 }, (|z 2 | > 1) 
z fo : {0,«f} {0,4}, (|z 2 |<l) {0,z 2 *}, (|za|>l) 

Table 1: The location of saddle points of q a , gt> are z a , Zt, correspondingly. The Z\ and z 2 are given in (|53p. 



The location of the solutions on the complex z-plane depends on the values of 5 and b: there are actually three 
distinct regions in the (b - 5) plane (see the left hand side of figure [T]) . We listed in table [1] the solutions of the 
saddle equations which correspond to the different values of (6, 5). In the table z\, z 2 denote 



Zl 



1 

Yb 



i + y/8bS - (2b - 1) 



Z2 



1 

2b 



1- V(2&-1) 



865 



(53) 



It is interesting that, the "trivial" z = saddle point is present in each region. However, when one continues 
analytically also the integrand of I 2 into the complex plane one finds that its contribution exactly cancels the 
contribution of the z — saddle point of I\. Hence it is enough to handle I\ with its non trivial saddle points. 
The integration contours which are displayed on the right hand side of figure [T] correspond to region I. 

For the actual evaluation with help of the saddle point technique let us recall that its general strategy tells us 
that we have to deform the integration contour of I\ in order the integration path passes over the saddle point 
(see the r.h.s. of figure [1} because the dominant contribution of the integral comes from the neighbourhood of 
the saddle point. Near the saddle point q a ^(z) can be approximated by their second order Taylor expansions 
around the saddle point, 



1 



q{z) = q(z s ) + -q | z< (z - z a )" 
where z s denotes the corresponding saddle point and 



q( z s) - d—, 



D = \q"(z s )\, z = z s +rexp 



-i / f lmq"(z s 

— arctan — ; — 

2 V VReg"(z s 



T 7r 



t 6 R. 



(54) 



(55) 



Practically speaking, the above form of q a ,b provides a parameterization (z{t)) of the integration path near 
the saddle point in which the argument of the exponential function have constant imaginary part along the 
parameterized integration contour (stationary phase). Using the above parameterization for an integral of type 
I\, one can write 



J dz^/zexp(nq(z)) « e M{Zs) J dr— ^/z(t) exp ( -fxD 



dr 



At exp I — jiD 



= \ z,e 



/ -2vr 



(56) 



where in the first line the notation J z means that the integration range is a small interval near z s , (r = 0). In 
the second line this range is extended to [—00, 00] because the contribution to a Gaussian integral is small far 
away from the saddle point when fj, is large. In addition we write ^pz~ s instead of yfz since it is a slowly varying 
function near the saddle point. Higher order terms of the Taylor expansion of *J~z give 0(1/ ji) corrections. 

It can be seen from (f57>| that the order of magnitude of the integral 1\ is determined by the real part of 
qa,b(z s ), which can be written for the regions i = (I, II, III) in the form q l a (z l s ) = K l + iT l , q l b (zl) = R l — iT l . 
After some algebra, where one pays attention to the branch cuts of exp(fj,q(z s )) appearing in (I56|) one can obtain 
the following explicit expressions for the real and imaginary parts: 



1 

4b 



6-1-25 26 + 1 



In (1 + S) +S In VP + S, 



ln(2) - -L (2b - 1 - y/(2b - l) 2 - 8W) + 



(5+l)ln(26+l- y/(2b-l) 2 ~ 8bd) -51n(26-l- (2b - l) 2 - 865), 



R 



in 



ln(2) - — [2b - 1 + v/(26- l) 2 -86,5 



+ (6 + 1) In (26 + 1 + y/(2b - l) 2 - 86(5) -5 In ( 26 - 1 + ^(26 - l) 2 - 865 ) , 



(57) 



(58) 



(59) 



9 



^/m - (26 - 1) 

46 



28 



latanV^-^- 1 ^ 



26 - 1 + 465 



latan^-^- 1 ) 2 
2 26 -1-26(5 



r/IJ 



0. 



(60) 
(61) 



Using these expressions together with the form of the normalization factor N from (|46[) and the exponential 
factor from (|4"T|) one obtains that 



< - 6/2 w O(l), therefore TV exp U(R 



ii, in 



6/2)) « 1, 



(62) 



hence J?^(/x/26) is exponentially small when 6, 5 are in regions II or III. Unlike in region I, where the exponential 
factors are cancelled, and the result is surprisingly simple: 



N' exp ( R 1 



1 1/4 



1 



(63) 



Using (f6"3")l in (|56[) together with the explicit form of q"(z s ), one can obtain the expression (j3"2"|) . This form of 
the asymptotic behaviour of the Laguerre polynomials is in agreement with the results of [16] , where asymptotic 
analysis of the Laguerre polynomials was derived with help of the generating function of these functions. 
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